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Maxwell's equations are solved for a periodic sequence of lens-like focusers 
separated by gaps. Each focuser consists of an arbitrarily thick slab of 
dielectric in which the dielectric constant tapers off radially with different 
quadratic laws in two perpendicular directions. Since there are no limita- 
tions on the thickness of the slabs, the solutions cover the complete gamut 
from a sequence of infinitely thin lenses with astigmatism to a continuous 
dielectric waveguide, and from spherical to cijlindrical lenses. 

The field configurations of the 77iodes and their propagation constants, as 
well as the transmission and cutoff bands, are calcidated. Any arbitrary 
input field distribution can then be expanded in terms of the normal modes, 
and the expansion determines the field everywhere. 

Foi-mulas derived for sequences of weak lenses turn out to give very good 
results even for lenses whose thickness and separation are equal to the focal 
length. 

I. INTRODUCTION 

One possible long distance transmission medium for optical waves 
consists of a periodic sequence of converging lenses. In order to negotiate 
unwanted but unavoidable bends of the axis of the sequence it is neces- 
sary to space the lenses as closely as possible. 1 Nevertheless, ordinary 
dielectric lenses exhibit substantial surface scattering, and therefore 
the minimum spacing between lenses depends on the tolerable transmis- 
sion loss. 

D. W. Berreman has shown that an effective lens can be made using 
gas with thermal gradients, 2 - 3 thus avoiding the solid-to-gas transition 
problems. 1). W. Berreman mid S. E. Miller 4 proposed a gaseous lens 
consisting of a tube with hot walls through which a mild gas current at 
lower temperature is forced to (low. At any cross section the tempera- 
ture increases from the center to the wall. The density and consequently 

2887 



2888 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1964 

the dielectric constant is then maximum on the axis and decreases 
radially roughly with a square law. Without the problem of scattering 
at the interfaces, tubular gas lenses can be closely spaced and the gaps 
may be comparable to the thickness of the lenses. 

The advent of such a new transmission medium makes it opportune 
and important to generalize the theory of modes in a sequence of thin 
lenses by determining the normal modes in an idealized structure which 
consists of a periodic sequence of arbitrarily thick slabs of dielectric 
whose dielectric constant tapers off radially with quadratic law. 

The preferential direction of gravity creates convection currents that 
may introduce astigmatism in the gaseous lenses. Such an aberration 
is included in our model by making the radial quadratic law of the 
dielectric different in two perpendicular directions. 

We calculate the modes of propagation of the idealized structure 
without including the solid walls surrounding the medium. Taking them 
into account would perturb the modes only slightly, introducing diffrac- 
tion losses. Just as in the case of a waveguide with perfect metallic 
walls, the idealized modes considered here are not attenuated, but their 
discussion is similarly expected to be useful in approximating: (a) the 
propagation constants: (b) the range of dimensions over which trans- 
mission is permitted or forbidden; (c) the extent of mode conversion at 
discontinuities or imperfections; and (d) the field at any point due to an 
arbitrary input such as an off -axis or tilted beam. Of these, (a) and (b) 
are treated in this article. 

The calculations are general enough that by changing the lens param- 
eters and the length of the gaps it is possible to cover uninterruptedly all 
the range from a sequence of thin lenses 5 ■ 6 ' 7 - 8 to a continuous dielectric 
guide, 1 ,J ' ,n and from spherical to cylindrical lenses. Up to now only the 
extreme cases, that is, thin lenses or dielectric guide and spherical or 
cylindrical lenses, have been considered in the literature; this article 
bridges the gaps. 

II. DESCRIPTION OF THE PROBLEM 

Consider a periodic sequence of dielectric slabs, shown in Fig. 1. 
The refractive index v of each slab is independent of z, but varies with 
different quadratic laws in the x and y directions as 

The refractive index n on the z axis and the characteristic parameters 
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Fig. 1 — Periodic sequence of arbitrarily thick and astigmatic lenses. 

of the lens-like medium, Li and La , permit adjustment of the parabolic 
distributions. The physical significance of L\ and L 2 will be treated 
below. 

In spite of the fact that it reaches negative values for large x or y, 
this dielectric distribution is useful because it matches the dielectric 
distribution of the gaseous lens, especially for small values of vx/L\ 
and iry/L-, . Besides, it turns out that the Held of most modes is negligi- 
ble in the region where the dielectric constant is small or negative, and 
consequently that region does not contribute essentially to the guidance 
of the modes. 

In the Appendix we solve approximately Maxwell's equations. The 
sequence of lens-like focusers supports hybrid modes EII p(/ , charac- 
terized by the indexes p and q. These integers indicate that the intensity 
of each transverse field component passes through p zeros in the x 
direction and q zeros in the y direction. 

The only approximation in the solution of Maxwell's equations 
consists in neglecting terms of the order of pX/Li and q\/L> compared 
to unity. X is the free-space wavelength. Typically X/Li and X/L 2 are 
of the order of 10~ 5 ; therefore, except for very high-order modes (p 
and/or q very large), the results must be satisfactorily precise. 

The modes have no electric field in the y direction nor magnetic field 
in the .r direction. The remaining components - E, , E. , 1I (/ and H_. 
in the dielectric, slabs and E z „ , Iv„ , TI M „ and H 2fl in the gaps - are found 
assuming as normal modes only those field configurations that repeat 
themselves periodically at each lens. Therefore the equiphase surfaces 
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of each mode are planes at z = and f = 0, as shown in Fig. 1. We 
reproduce here only E x (55) and E x „ (58). 

All the other components can be deduced from them with the help 
of (45). 

E, = e*p {-, [„. (, - 4 - Q - (, + I) tan"' g ta g 

-hms-e)]-©'-©) « 

■H,(v^2)H.(vSj)i 
•H,(V2 ±) H, (V2 i) 

\ Pl/I/ \ Po2/ 

where [see (55) to (67)], fc = w\/eoM = 2tt/X is the free-space prop- 
agation constant and H M (a) is the Hermite polynomial of order /*• 

The physical significance of the symbols Si , s 2 , s ffX , s„ 2 , Ri , R 2 , etc. 
will be developed below. We give first their mathematical meaning and 
in order to avoid repetition, from now on the letter m will stand for 
either the subindex 1 or 2, depending on whether the symbol under 
consideration refers to a dimension in the plane y = or x = respec- 
tively. Calling the thickness of each dielectric slab /, and the gap between 
them b, 

s „ = „. ( ; + g- "■» *■■■ )' (4 ) 

\1 — C„, tan <p„J 

Sy m = W m {\ + C m CtnpM)*(l — C m tell Ipm)* (5) 

(6) 
(7) 
(8) 
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-•JFeR-^T-S (l0) 

>-=«»V i+ (*!?)' (12) 

Let us find the physical significance of i? m , R gm , p m , p om , s m , s gm , 
w m and L m . Equating in (2) and (3) the imaginary parts of the expo- 
nents to constants we obtain two equations of equiphase surfaces (wave- 
fronts), one applicable within a lens and the other in a gap. At the z 
axis, each wavefront has a radius of curvature in the plane y = which 
in general is different from that in the plane x = 0. Within a lens those 
main radii of curvature are Ri and R 2 [see (9)], while those in a gap are 
R i and R g2 [see (11)]. If L x = L 2 , then R t = R 2 and R gi = R g2 . 

For the fundamental mode p = q = 0, at a given abscissa z or f the 
field amplitudes (2) and (3) decrease with different Gaussian laws in 
the x and y directions. The distances at which the field is 1/e of the 
maximum occurring on the z axis are the beam sizes pi and p 2 [see (10)] 
within a lens, and p g i and p g2 [see (12)] in a gap. 

For 2 = and £ = we find from (10) and (12) that p m = s m and 
P„m = s gm . Therefore s,„ and s gm are the beam sizes at the planes of 
symmetry of each lens and each gap respectively. 

The physical significance of w m becomes obvious on reducing the gaps 
between lenses to zero. Then instead of a sequence of lenses we have an 
uninterrupted dielectric waveguide and we derive from (7), (4), (5) 
(10) and (12) that 

Pm = Pum = S m = 8 gm = W,„ . (13) 

Therefore in the continuous guide the propagating normal modes do 
not change size along z, and for the fundamental mode Wi and w 2 measure 
the beam sizes in the x and y directions. 

From (10) we find that within a lens the beam sizes p\ and p 2 in the 
y = and x = planes vary periodically along z; their periods are L\ 
and L 2 respectively. 

For the particular case in which Li = L 2 the field in the gap (3) 
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coincides with that found by Boyd and Gordon 6 for the resonator made 
with confocal mirrors of infinite aperture. 

III. TRANSMISSION AND CUTOFF CONDITIONS 

Both s m [see (4)] and s gm [see (5)J must be real quantities, otherwise 
the fields given in (2) and (.3) become infinite as x or y — * oo. This 
establishes that a mode can propagate in the sequence of lenses either 
when 

C m ^ eta*. (14) 

or when 

C m ^ -tan*™. (15) 

Their equivalents in explicit form are 

b ^ (2L m /nir) ctn (irt/2L m ) (16) 

and 

b ^ -(2L m /mr) tan (irt/2L m ). (17) 

Which equation must we use? Since b and L„, are positive, (14) or (16) 
must be used when <p m = vt/2L m falls in an odd quadrant and (15) or 
(17) when it falls in an even quadrant. Naturally, if these equations are 
satisfied for only one of the two indexes, that sequence of lenses cannot 
propagate any nonattenuating mode. 

If b = 0, the sequence of lenses is reduced to a continuous waveguide 
and transmission takes place, as it must, no matter what the values 
of <pi and <p 2 are. If now we increase the gap 6, transmission will take 
place as long as (16) or (17) is satisfied. 

IV. DISCUSSION OF THE FIELD INSIDE AND OUTSIDE THE LENSES 

The sequence of lenses admits a complete set of modes. For each mode, 
the field inside (2) and outside (3) the lenses is a wave traveling in the 
z direction whose amplitude, period and equiphase surfaces (wave- 
fronts) vary along z. 

The amplitude depends on .r as a product of a Gaussian function and 
a Hermite polynomial (parabolic cylinder function) whose degree 
depends on the mode under consideration. A similar type of variation 
occurs along y. 

In Fig. 2 we plot qualitatively the beam sizes p„, and p gm for ip m = 
irt/2L m in the first, second and third quadrants. For <p m in an odd quad- 
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rant, as in Figs. 2(a) and 2(c), the maximum and minimum beam sizes 
within each lens are 



= s m = w. 



and 

2 



w, 

— - = w, 

Sm 



(\ + C m Ctll ip m \ / 1q n 

{l-Cmtan^) U8) 

/ l-^tan^ V ( } 



The period between two successive maxima is L m . The square root of 
the product of the maximum and minimum beam sizes in the dielectric 
is a constant 

(pm max Pm min) = W m 

and coincides with the beam size w„, of the lens-like medium. 

In the gap, the only extremum for the beam size is a single minimum 
which occurs at f = and, from (5) and (12), corresponds to 

Pom min = S gm = W m (l + C„, Ctll *.)* (1 - C m tail <p m )\ (20) 

If <p„, = irt/2L m falls in an even quadrant, as in Fig. 2(b), the mini- 
mum and maximum beam sizes interchanged from the odd quadrant 
are (18) and (19) respectively. Again (20) corresponds to the unique 
minimum in each gap. 

V. SPECIAL CASES 

Let us consider the field in a gap assuming 

L, = U = L 

and 

t/L = v or tp x = ipi = t]{tt/2) (21) 

where rj is an integer. Then unless the gap b = 0, the minimum beam 
size in the gap p gni min (20) becomes infinitely large and the electric 
field (3) is reduced to a plane wave travelling in the z direction. If more 
generally only <p\ = r)(ir/2), but <p 2 is unrestricted, then the wave fronts 
are cylindrical surfaces parallel to the .r axis. 
Consider again 

U = U = L 

but 




Fig. 2 — (a) Beam size for <p m = - ■=— in first quadrant, t/L m < 1; (b) beam 

& Li m 

size for <p m = - j— in second quadrant, 1 < t/L m < 2; (c) beam size for <p m = ■= j— 

in third quadrant, 2 < t/L m < 3. 
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& = C 2 = ctn<pi = ctn<p 2 (22) 

or 

d = d = — tan ^1 = — tan ¥> 2 • (23) 

Then according to (20) the minimum beam size p„i mi „ = p 2 min = 
and the field in the gap (3), for p = q = 0, becomes 

E Ia = exp -ikt (l + ^~) ■ W 

The wavefronts close to the f axis are concentric spheres and their 
centers coincide with the point x = y = f = 0. 

Therefore the two conditions indicated above correspond either to 
plane waves in the gap or to concentric waves (if one observes only the 
field in the region close to the f axis). They are equivalent to those in 
Fabry-Porot resonators with plane and concentric mirrors. ' 

The condition under which the beam is closely concentrated on the 
z axis is found by minimizing the maximum beam size within a lens, 
s„, [see (18)] or wj/s m [see (19)] depending on whether ip m is in an odd 
or even quadrant. 

If the gap 6 decreases, the value of s m or w m 2 /s m also decreases; for 
6 = the sequence of lenses becomes a dielectric waveguide, the beam 
size does not vary with z, and its value is s m = wj/s m = w m . On the 
other hand, if the thickness t of each lens is the only variable, the mini- 
mum of s,„ or w m /8 m is achieved by making 



oi- 



— - = — - = if ip m is an odd quadrant 

dt bipm 



1 d 1 n ., . 

— = if <p m is m an even quadrant. 



(25) 



dt S,„ d<p„, S„ 
These conditions lead to the same requirement, namely: 

C m = ctn2v>m (26) 

or its equivalent 

b = (2L m /nir) ctn (irt/L m ) (27) 

which, replaced in (18) or (19), determines the minimized value of the 
maximum beam size within each lens 

S m ,nin = (Wj/Sm) min = W m [(l + Cj)* + Cj. (28) 
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For the same condition (26) or (27), the beam size in the gap at any 
abscissa f is derived from (5), (12) and (26) 

P„„, = W .(l + 0'[ 1+r ^(f)]'. (29) 

VI. SEQUENCE OF WEAK ASTIGMATIC LENSES 

Before considering weak lenses, let us relate the characteristic lengths, 
Li and L 2 of the lens-like focusers to their focal lengths in the planes 
y = and X = 0. To calculate the focal length in the y = plane (see 
Fig. 3) the ray trajectory is determined from the equation 

d'x/dz = (1/p) {dv/dx). 



Taking the refractive index v from (1) 



(fx 



1 



^ Vl - {irx/UY dx 
For paraxial rays 

tx/Li « 1 



^Vi - (jx/uy. 



(30) 
(31) 

(32) 



and within a lens the trajectory of a ray entering parallel to the z axis 
at a distance x is 

X = .To cos (wz/Li). 

The angle of refraction at the output surface is 

d T = (mr/L^Xo sin (irt/Li). (33) 




Fig. 3 — Ray trajectory in the plane y = of a lens. 
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Then from simple geometric considerations in Fig. 3, the focal length 
/i results 

f x = h. . (34) 

ten sin (wl/Li) 

We assume now weak lenses. They are characterized by 

<p,„ = (tt/2) (f/L,„) « 1, (35) 

and in all previous results each circular function can be replaced by its 
leading term. 

Because of the inequality (35) the characteristic length of the focusing 
medium L, in (34) can be calculated explicitly by 

U = irVWi • ( :5(i ) 

Similarly, for the plane x = 

I, = rVWt . (37) 

The weak lens requirement (35) then becomes 

Using (36) and (37) together with the simplifying assumption (38) we 
re-evaluate the maximum and minimum beam sizes (18), (20) for 
weak lenses dp,,, in first quadrant), 

'*/ m \ 2 \7 1 + (nb/t) V 



s,„ = 



Mr- J \\ - (b/4f m ) 



(39) 



The distance h between the principal planes may be of interest. Using 
(33) and (34), this distance turns out to be 

hm =' 2 Lm t an^-(. (41) 

Til ZlJm 

Expanding the circular function in series, keeping only the first two 
terms and substituting L,„ by their equivalents (30) and (37) we ob- 



. • 14 

tarn, 



h ' I " ~ * ) + T^T ■ (42) 
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6.1 Example 

Let us assume a sequence of gaseous lenses such that 

b = t = ji = h = 0.25 m 

X = 0.6328 10~ 6 m 

n « 1. 

For these dimensions v?i = <p 2 = 0.5, and therefore the weak lens in- 
equality (38) is hardly satisfied. Nevertheless, let us go ahead and 
calculate extreme beam sizes s x = s 2 and s g i = s r;2 as well as the charac- 
teristic length L\ = L 2 of the lens using (39), (40) and (36) 

Sl = s 2 = 0.286 mm 

s„i = v = 0.248 mm (43) 

Li = L 2 = 0.785 m. 

Let us calculate again the extreme beam sizes using the exact expres- 
sions (18), (20), deriving L from (34) 

si = s 2 = 0.276 mm 

s g i = s„., = 0.224 mm (44) 

L, = L, = 0.704 m. 

The two sets of results (43) and (44) are reasonably similar and show 
the usefulness of weak lens formulas even for lenses with comparable 
values of t and /. 

VII. CONCLUSIONS 

The properties of the modes in a sequence of thick, astigmatic and 
unbounded lens-like focusers are similar to those in a sequence of thin 
infinitely large lenses. 

The modes are hybrid and described by parabolic cylinder functions 
(product of Gaussians times Hermite polynomials). Transmission takes 
place as long as the gap between lenses is smaller than a value given in 
(16) or (17). 

The maximum beam size can be reduced by decreasing the distance 
between dielectric slabs. Nevertheless, if the gap is fixed, the minimiza- 
tion of the maximum beam size can be obtained by selecting the dielec- 
tric properties or the thickness of each focuser according to (27). 

Simplified formulas derived for sequences of weak lenses yield good 
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approximations oven for lenses whose thickness, separation and focal 
length are comparable. 
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APPENDIX 

Solution of Maxwell's Equations in a Sequence of Thick Astigmatic Lenses 

We will obtain, first, a general enough solution of Maxwell's equations 
for one of the lenses; see Fig. 1. Then by making n = 1 and Li = L 2 — » » , 
we will deduce a general solution for the uniform gap between lenses, 
and finally we will match the tangential fields to satisfy the boundary 
conditions. For modes with only four field components, E x , E z , H„ and 
K z , Maxwell's equations become 

dE x dE,_ . u 

OZ ox 



dE x 

di! z _ dilj 

dy c)z 

dx 
whore /x, the magnetic permeability, is a constant; 



= — /weEj- 
= — jueE z 



(45) 



— = n 



1 - fey - ©i ^ 

e, L, and L, are arbitrary constants; and w-y/eoM = 2x/X = k is the 
free-space propagation constant. 

By eliminating variables and by neglecting terms of the order of 
\/Li and \/L* as compared to unity we obtain identical equations for 
E, and H„ . For E, , 



d 2 E x d 2 E x a 
dx 2 dy- dz* 



% + ««r[i -($-(%)']*.-<>. w 



* In practice X/Li and \/L» are of the order of 10 6 . 
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This equation is separable and a general solution is 



X 2^2^ A VII exp -iknz ( 1 -!— ^— — ) 

„=o m-0 L \ Zn Li Zn Li/ . 

\ W\ J \ Wi ) 



(48) 



where v and ^ are integers, and A vll is an arbitrary constant. Using m 
to indicate either subscript 1 or 2, 



7T 



(49) 



The function 



HM) = (-i)V WW* 



is the Hermite polynomial 12 of order v. Hermite polynomials of lowest 
degree are #„(£) = 1; //,(£) = 2£; ff 2 (£) =4^-2; and ff,($) = 
8f 3 - 12$. 

Expression (48) can be simplified provided that the important terms 
of the summation are those for which 



and 



vk/Lx « 1 



n\/L 2 « 1. 



(50) 



Then the square root in the exponent can be replaced by the first two 
terms of a power series expansion and 



/ ? P ( - —. i - —J |E-1, exp [t(WLi)]W, 



»'i- 



«'•> 



(£')] 



X 



[g 5„ exp [*(WIa)]#, (^)] 



(51) 



where 



= hnz 



1 - 



4« 



fe+ai- 



We will look for a periodic field configuration that reproduces itself 
at each lens. For reasons of symmetry, then, the planes of symmetry of 
the lenses (z = 0) and gaps (f = 0) must be equiphase surfaces. 



MODES IN LENS-LIKE FOCUS ERS 

We choose the field at the plane z = to be 



2901 



< I [ — '— , — "— ) H, 



V2A H (V2V 



syy " \ Si 

where p and q are integers and si and s 2 are arbitrary parameters for 
the time being. Therefore for z = 0, we obtain from (51) 



exp 



-5 



-—AH 



and 



£a. 



^^[-h) H 



V2a 



V^2/ 



exp 



exp 



*r 



//, 



(52) 



(V2u 



.,H.\^l (53) 



Using the orthogonality properties of the Ganssian-Hermitian product, 1 
we obtain 



a,- v? r exp \-eU + ±X\ 



•*>(^W^)« 



(54) 



and a similar expression for J5„ . Replacing the result in (51) and per- 
forming, as in Ref. 13, first the summation in v and /x and then the 
integration in £, the transverse field component inside a lens expressed 
in closed form results 

E, = exp j-i [fcn [z - JL - ^Q - (v + ^) 

-- i g-s)-G + S-- , & , -H)] 
x -[-©'-©> K)*«K) 

where, for m = 1 or 2 



(55) 



R m = — 

IT 



1+ — 



fey 



. io„, \ \ . iswz 

1 'id 8in ZT 



1- ctn — 



(56) 



and 



2902 THE BELL SYSTEM TECHNICAL JOUENAL, NOVEMBER 10G4 



P,„ = s„ 



The electric field in the uniform dielectric gap between two lenses 
can be derived from the previous expression by making 

n = 1 

and 

1/L m = 

and by substituting another symbol, s om , for s m . Again we demand the 
plane of symmetry of the gap, z = (b + t)/2 (see Fig. 1), to be an 
equiphase surface. This is achieved by substituting f = z — (b + t)/2 
for 2. 

The electric field in the gap is then 

-R) --a-©-©} « 

•H P (V2-W(V2^-) 

\ Pol/ \ Pal/ 



where 



and 



th am 



Pg» 



¥b+(£f] 



1 + (sbT- (fi0) 



To match the fields (55) and (58) at the interfaces, the x and y 
dependences of the field at both sides must coincide. The fact that it 
can be matched guarantees that Maxwell's equations are satisfied 
simultaneously in lenses and gaps. It can be verified that if the tangential 
electric field continuity is satisfied, the tangential magnetic field con- 
tinuity is also guaranteed. By considering waves propagating in both 
directions, it could be possible to take into account reflections at the 
interfaces, but we shall instead assume that at each interface there is a 
matching mechanism that prevents reflections. Notice that in the case of 
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gaseous lenses the small changes of dielectric constants automatically 
insure negligible reflection at the interfaces. 

The exact matching of the fields at the interfaces is achieved by 
making equal the coefficients of x, ij, .r 2 and if in both expressions (55) 
and (58) at the boundary z = 1/2 of the lens and f = -6/2 of the gap. 
Then 

Rm{t = t/2) = Ro<»<.S --6/2) ("1) 

Pm(t - »/2) = PgmlS --(6/2)1 ■ ( 62 ) 

From them, together with (56), (57), (59) and (60), we deduce the 
values of s m and s um that guarantee the matching at the interfaces. 
They are : 



= w. 



|~ 1 + C m ctn y m ~l* ^ 

Ll — C m tan tp m \ 



and 

s gm = w m (l + ('„, ctn v> m ) { (1 - C m tan ^ m ) ! (64) 



where 



("„, = n(*/2) (6/L m ) (65) 

^„, = (tt/2) (£/L.) (66) 

u) w = (1/tt) VxIZM- (67) 
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